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Convective storms can pose a serious risk to aviation operations since they are often 
accompanied by turbulence, heavy rain, hail, icing, lightning, strong winds, and poor 
visibility. They can cause major delays in air traffic due to the re-routing of flights, and by 
disrupting operations at the airports in the vicinity of the storm system. In this study, the 
Terminal Area Simulation System is used to simulate five different convective events 
ranging from a mesoscale convective system to isolated storms. The occurrence of 
convection induced turbulence is analyzed from these simulations. The validation of model 
results with the radar data and other observations is reported and an aircraft-centric 
turbulence hazard metric calculated for each case is discussed. The turbulence analysis 
showed that large pockets of significant turbulence hazard can be found in regions of low 
radar reflectivity. Moderate and severe turbulence was often found in building cumulus 
turrets and overshooting tops. 
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Nomenclature 

decibels of radar reflectivity, Z 
Mean Sea Level 
Next Generation Radar 
acceleration due to earth’s gravity 
atmospheric pressure 
vertical component of velocity 
atmospheric temperature 
Cartesian coordinate in east direction 
Cartesian coordinate in north direction 

Cartesian coordinate in vertical direction, altitude above ground 

aircraft’s normal load acceleration 

grid size along x-coordinate 

grid size along y-coordinate 

grid size along z-coordinate 

averaging interval along x-direction 

averaging interval along y-direction 

variance of vertical velocity 

root mean square of aircraft’s normal load acceleration 


I. Introduction 

A ircraft encounters with turbulence are the leading cause of in-flight injuries (Tyrvanas 2003), and on a few 
occasions have resulted in passenger and crew fatalities. According to a National Transportation and Safety 
Board (NTSB) report posted on the internet, there were 19,562 aircraft accidents between 1994 and 2003 and 
weather was a factor in approximately 21% of these accidents. In the Federal Aviation Regulations (FAR) Part 121 
(Air Carrier Operations) category, seventy four percent of the weather related accidents were due to turbulence 
encounters (http://www.asias.faa.gov/aviation_studies/weather_study/studyindex.html). A substantial number of 
these turbulence encounters was likely due to convection induced turbulence (CIT). Studies have shown that most 
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of the injuries were caused by sudden and unexpected encounters with turbulence that occurred within or in 
proximity to convective activity (Kaplan 2005). Encounters with severe CIT are usually brief, lasting several 
seconds to a few minutes. Severe encounters may occur as aircraft attempt to navigate around high radar reflectivity 
regions or when cells rise from below the flight path, especially when visibility is poor. 

In the fall of 2000 and the spring of 2002, two flight campaigns were conducted by the National Aeronautics and 
Space Administration (NASA) to test the feasibility of a radar-based airborne turbulence detection system. NASA 
Langley Research Center’s B-757 research aircraft was flown into regions of convectively induced turbulence. 
Turbulence encounters were typically less than 30s in duration and coincided with the time spent within the cloud 
(Hamilton and Proctor 2002; Hamilton and Proctor 2006). Furthermore, the events occurred outside of moderate 
and high reflectivity regions (> 35 dBZ) which were avoided during the flight experiments. Strong turbulence 
encounters occurred with radar reflectivity factors as low as 0 dBZ and no correlation was apparent between 
turbulence intensity and the radar reflectivity factor (RRF). 

In commercial aircraft encounters with CIT, onboard weather radar may fail to detect regions of low RRF, and 
the incidents are therefore classified as encounters with clear air turbulence (CAT). For example, on 15 July 2010, a 
Delta Airlines’ Boeing 767-300 encountered turbulence 250 nmi north-northwest of Guam while at 36,000 ft. The 
onboard weather radar indicated nothing prior to the incident. During the encounter, the aircraft underwent vertical 
accelerations between +1.5g and -0.3g for approximately 5s. A flight attendant walking down the aisle fell down 
and fractured her ankle. The investigation into the incident concluded that the probable cause of the accident was an 
inadvertent encounter with clear air turbulence (NTSB Identification: DCA10FA076). This conclusion is suspect 
since the crew reported scattered cloud build ups around the aircraft at the time of the incident, and synoptic weather 
conditions did not support the likelihood of CAT. 

Simulations with the Terminal Area Simulation System (TASS) have been used in the past for the development 
and testing of algorithms used by airborne turbulence detection systems and in defining turbulence hazard in terms 
of aircraft-centric metrics (Hamilton and Proctor 2006). These simulations also provide a detailed insight into the 
formation and evolution of convective systems. In the following sections, the TASS, the Weather Research and 
Forecast (WRF) mesoscale model, and the turbulence hazard metric used in this study are briefly described. The 
paper examines five severe convective events simulated with TASS. The simulation results are compared with the 
Next Generation Radar (NEXRAD) data (http://www.roc.noaa.gov/WSR88D/) and other observations. Aircraft- 
centric turbulence characteristics associated with each convective event are also discussed. 

II. Terminal Area Simulation System (TASS) 

TASS is a state-of-the-art, cloud-resolving, large eddy simulation atmospheric flow model (Proctor 1987, 1996, 
2002; Proctor et al. 2002). TASS computes the primitive variables non-hydro static equations in three dimensions 
and is capable of resolving flows at multiple spatial (grid resolutions varying from less than \m to 2 km in the 
horizontal) and temporal scales (few seconds in the case of turbulence eddies to hours for long-lived convective 
phenomena). The model solves prognostic equations for potential temperature, water vapor, cloud droplets, ice 
crystals, rain, snow and hail, and includes a microphysics package for cloud and precipitation development (Proctor 
1987). Subgrid scale diffusion is parameterized via a Smagorinsky-type turbulence closure (Smagorinsky 1963) and 
surface layer processes are computed based on the Monin-Obukhov similarity theory (Stull 1997). The TASS model 
equations are discretized using quadratic-conservative fourth-order finite-differences in space for the calculation of 
momentum and pressure fields (Proctor 1996) and the third-order Leonard scheme (Leonard 1995) is used to 
calculate the transport of potential temperature and water vapor. A Monotone Upstream-centered Scheme for 
Conservation Laws (MUSCL)-type scheme after van Leer (van Leer 1979; Ahmad and Proctor 2011) has also been 
implemented in TASS for the transport of water substance variables. The TASS computational mesh uses the 
Arakawa C-grid staggering (Arakawa and Lamb 1977) for specifying velocities and thermodynamic quantities. The 
Klemp-Wilhelmson time-splitting scheme (Klemp and Wilhelmson 1978) is used for computational efficiency in 
which the higher-frequency terms are integrated by enforcing the CFL criteria to take into account sound wave 
propagation due to compressibility effects. The remaining terms are integrated using a larger time step that is 
appropriate for anelastic and incompressible flows. Nonreflecting Orlanski boundary conditions (Orlanski 1976) are 
imposed on open/outflow lateral boundaries. Options also exist for periodic lateral boundary conditions. A sixth- 
order filter is used to damp-out spurious oscillations in the velocity field that may arise due to the use of centered- 
differencing of momentum and pressure terms. Filter and sponge conditions (Perkey and Kreitzberg 1976) are 
applied to the top four levels to minimize flow interactions and gravity wave reflections with the closed top 
boundary. The model code is written in FORTRAN and has been fully parallelized for distributed computing 
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platforms using the Message Passing Interface (MPI). Excellent scalability and vectorization has been demonstrated 
on high-end supercomputing platforms such as NASA’s SGI Altix cluster. 

III. Initialization of TASS Simulations 

Two steps were needed for achieving high resolution simulations of each of the five cases described in this 
paper. In the first step, a numerical weather prediction model was run in order to obtain representative 
environmental conditions. In the second step, the environmental conditions were used as initial conditions in TASS 
to simulate the detailed structure and evolution of the convective event. The environmental conditions were 
obtained with the National Center for Atmospheric Research’s (NCAR) Weather Research and Forecast (WRF) 
model (Klemp et al. 2000). The WRF model has been developed by NCAR with help from various partners in 
academia and government. WRF is a mesoscale weather prediction model with governing equations for the time- 
dependent, three-dimensional, nonhydro static, fully compressible Navier-Stokes equations. In the vertical, WRF 
uses a terrain-following pressure formulation, where the top of the domain is a constant pressure surface. The model 
grid structure is based on Arakawa C-grid staggering (Arakawa and Famb 1977). Both horizontal and vertical 
advection terms can be discretized anywhere from second- to sixth-order spatial accuracy. Fifth-order upwind 
differencing in the horizontal direction and third-order upwind differencing in the vertical direction are 
recommended by the developers of the WRF model. The time integration is done using an explicit third-order 
Runge-Kutta time marching scheme, with smaller time steps for acoustic and gravity wave modes. Extensive 
physics packages for modeling the boundary layer turbulence, cumulus parameterizations, cloud microphysics, 
radiation physics and surface layer interactions are included in the model distribution. Fligh resolution terrain and 
other land use datasets are part of the standard WRF distribution. Numerical forecasts were conducted with WRF to 
provide initial soundings for each of the five cases described in this paper. The initial environment in TASS is 
assumed horizontally homogeneous, and varies only in the vertical condition in accordance with the sounding 
extracted from the WRF simulation. Convection with cloud and precipitation development in TASS is triggered by 
imposing an artificial bubble of warm air. Three-dimensional velocity, thermodynamic, and moisture fields evolve 
in the TASS simulation as the convection interacts with the environmental fields. 

IV. Flight Turbulence Hazard Metric 

The preferred aircraft turbulence hazard metric (Bowles and Buck 2009; Hamilton and Proctor 2006; Proctor and 
Hamilton 2005) is a five-second root mean square of the aircraft’s normal load accelerations (cr An ). This metric is 
preferred because: 

1) It is airplane centric 

2) The metric is easy to calculate 

3) It statistically quantifies the sharp bumps and accelerations that passengers feel when flying in an aircraft 

4) And it is understood by the airplane operators 

For a given turbulence event, the RMS of normal load accelerations (cr An ) is strongly correlated with the peak 
normal load (|A«|). This correlation is shown in Figure 1. Since the relationship of the peak normal load (|A«|) with 
the root mean square of normal load accelerations is known; the hazard levels can be assigned. The turbulence 
metric thresholds are subjective, but for warning pilots and dispatchers they can be categorized as follows (Bowles 
and Buck 2009): 


jap 

© 

VI 

R 

b* 

Light Turbulence 
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Moderate Turbulence 

0.2 < cr A „ < 0.3g 

Moderate to Severe Turbulence 

0.3 <a An < 0.6 g 

Severe Turbulence 

0.6 <<t a „ 

Extreme Turbulence 


The NTSB accident reconstructions shown in Figure 1 have values greater than 0.35g root mean square normal 
load acceleration (cr An ) or lg peak normal load (|A«|). Since aircraft response is most affected by the along-track 
gradients in vertical velocity, the aircraft normal loads can be estimated from aircraft flight dynamics models. A 
more general algorithm that can easily be applied to large data sets has been used by Proctor et al. (2002a, 2002b) 
and Hamilton and Proctor (2003). This method relates aircraft RMS normal load acceleration to the variance in 
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vertical velocity computed with a 1 km x l km moving box. The cr w fields can be computed for any horizontal plane 
by using a moving average: 


c w (x,y) 
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where, L x and L y , are the averaging intervals along the x and y coordinates respectively. The average vertical wind 
w , is computed from the vertical wind, w, as follows: 
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Figure 1. Relationship between the peak normal load (Aw) and the root mean square of normal load 
acceleration ((T An ) assuming a 5 second window. The dataset is constructed from measurements for 606 
turbulence encounters (Bowles and Buck 2009). 


The value for the averaging interval, L x = L y = lkm, is chosen which roughly corresponds to a five second 
averaging period for a commercial aircraft flying at cruise speeds. Turbulence intensities, cr An are estimated from cr w 
using the following relationship developed by Bowles and Buck (2009): 
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where, V is the aircraft’s true airspeed and W is the aircraft weight in klbs. The turbulence length scale, / is set to 
500m. The accuracy of Eq. (3) has been verified from NASA B-757’s in situ flight data. The parameters a(z ) and 
b(z) are designed for different aircraft types and calculated from look-up tables. a(z) is defined by the following 
equation: 


a(z) = c x + c 2 z + c 3 z 2 + c 4 z 3 + c 5 z 4 + c 6 z 5 + c 7 z 6 (4) 

where, for a Boeing 757, the constants, c i, c 2 , c 3 , c 4 , c 5 , c 6 , and c 7 are as follows (Bowles and Buck 2009): 
c, = 0.1367s _1 

c 2 =-4.759443 xl0“ 5 m _1 5 _1 
c 3 =3.453818xl0“ 8 m“V 1 
c 4 =-1.031272xl0 _11 m“ 3 5 _1 
c 5 =1. 457209 x10“ 15 ot“ 4 5 _1 
c 6 = -9.709422 x 10“ 22 m _5 5 _1 
c 7 =2.448279 xlO- 24 ^-^- 1 

The parameter b(z) is defined by the following equation: 

b(z) = c x + + c 3 z 2 + c 4 z 3 + c 5 z 4 + c 6 z 5 (5) 

where, for a Boeing 757 the constants, c u c 2 , c 3 , c 4 , c 5 , and c 6 are as follows (Bowles and Buck 2009): 

q = -0.0179545s -1 
c 2 =3.996083 xl0" 5 m _1 5 _1 
c 3 =-1.243271 xl0~ s m~ 2 s~' 
c 4 =1.760674 xlO^m-V 1 
c 5 = -1.153713x 10 _16 m _4 5 _1 
c 6 =2.807013 xl0 _21 m _5 5 _1 

The parameter F„(z) is calculated as follows: 

232.21 ms" 1 

2 

c 1 +C 2 Z + C 3 Z 

where, for a Boeing 757 the constants, c h c 2 , and c 3 are as follows (Bowles and Buck 2009): 

c, = 150.1493 ms~' 
c 2 =6.2023 x 10 _3 5 _1 
c 3 =3.0298 xl0 _7 m _1 s _1 

The algorithms described above are valid for mesh resolution on the order of 100m or less. Simulations of large 
convective systems at such high mesh resolutions require significant computing resources and very efficient codes. 
The simulations presented in this paper were run with mesh resolution of 100m or less in order to take into account 
the turbulence scales that affect the aircraft response (Proctor et al. 2002). 



if z > 9 A km 
if 3 km < z < 9 A km 
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V. Results 


In this section, the results of five TASS simulations are discussed in detail. The events ranged from a large 
mesoscale convective system, to a supercell hailstorm and small isolated convective cells (Table 1). The simulation 
results are compared with NEXRAD data and other surface observations. Aircraft-centric turbulence characteristics 
for each case are also presented. 


Table 1: TASS Simulations 


Case 

Event Type 

Location 

Ax (m) 

A 

Mesoscale Convective System 

Minnesota 

100 

B 

Single Cell Tropical Maritime 

Southwest of Puerto Rico 

60 

C 

Supercell Midwestern 

Indiana 

100 

D 

Multicell Tropical Maritime 

Northwestern Puerto Rico 

80 

E 

Multicell Midwestern 

Illinois 

80 


A. Mesoscale Convective System (17 July 2008) 

A Mesoscale Convective System (MCS) is an organized collection of several thunderstorms that can maintain its 
strength for several hours (Ashley 2003). There is no restriction on the shape of an MCS which can be elongated 
like a squall line or round/oval shaped like a Mesoscale Convective Complex (MCC). Maddox (1980) formally 
defined an MCC as an MCS with the following characteristics based on infrared satellite images (definition taken 
from the AMS glossary): 

• Size: Area of cloud top defined by an infrared derived temperature of -32°C or less: 100,000&m 2 or more 
(slightly smaller than the state of Ohio), and area of cloud top with temperature of -52°C or less: 50,000 km 2 
or more. 

• Duration: Size criteria must be met for at least 6 hours - Peak storm intensity normally at night. 

• Eccentricity: Minor/major axis at least 0.7. 

• Any large system of organized thunderstorms which does not meet the strict criteria of Maddox is generally 
characterized as an MCS. 


Due to their extensive size, MCS/MCC’s are difficult to simulate with an explicit cloud model. Past studies have 
examined MCS/MCC’s and their environment with mesoscale models (e.g., Zhang and Fritsch 1988). These studies 
usually assume grid sizes of several kilometers or more; and therefore, cannot explicitly resolve the MCS/MCC’s 
convective updrafts. The objective of this case was to simulate the most active parts of an MCS using TASS. The 
challenge was to maintain a sufficiently fine grid resolution in order to resolve the cumulus thermals and other 
important scales of motion within a computational domain large enough to encompass most of the MCS/MCC. 

Observations indicated that on 17 July 2008, a large system of thunderstorms moved eastward from Nebraska to 
southern Minnesota, northern Iowa and Wisconsin. The infrared imagery from GOES-E shows consistent 
temperatures below -32°C over a large area (Figure 2) indicative of intense convection. The system moved 
generally eastward at approximately \lm/s. The infrared imagery from GOES-E satellite shows that this system 
persisted for more than eight hours. Figure 3 shows the same system at 1500UTC with a clearly defined outflow 
boundary which is due to high pressure near the surface at the center of the MCS. The system produced large 
amounts of rainfall and hail (25 mm diameter hail observed near Wasesa between 1430UTC-1500UTC and 18mm 
hail observed at Blue Earth around 1536UTC) over southern Minnesota. 

Wind gusts as high as 22 m/s (50 mph) were measured at Freeborn, Minnesota and an EFO-tomado was reported 
at Amboy around 1655UTC. The structure of the low-level radar reflectivity field of the MCS was characterized by 
1) a comma shaped line of strong convective cells extending southwest along the southern boundary of the system 
(approaching the Iowa-Minnesota boarder); 2) small clusters of convective cells located along the northern 
periphery of the system; 3) a line of convective cells extending southeastward and eastward within the southeastern 
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portion of the system; and 4) a broad area of light rain covering the central region of the system (see Figure 4). 
Satellite photographs indicated that the MCS had an extensive, eastward spreading anvil (Figures 2 and 3). 



Figure 2. Surface weather map for the continental United States overlaid with GOES-E IR imagery at 
1200UTC on 17 July 2008. Colder temperatures indicate higher cloud tops. Source: National Centers for 
Environmental Prediction (NCEP). 



Figure 3. Surface weather map for the continental United States overlaid with GOES-E IR imagery at 10:00 
AM CST (1500UTC) on 17 July 2008. A zoomed view of the area covering Minnesota is shown in the inset. 
Source: National Centers for Environmental Prediction (NCEP). 
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1. Simulation Setup 

The TASS computational domain was set to 100 km x 100 km in the horizontal with a vertical depth of \lkm. 
Figure 4 shows the TASS model domain relative to the storm system over Minnesota. Since a long simulation time 
over a large domain was needed, the simulation was divided into two phases. In the first phase, a moderately coarse 
mesh with a grid size of 200m in all coordinate directions was used and the simulation was carried out until an MCS 
became established. After two hours and ten minutes of simulation time, the simulation was re-started with a re- 
gridded mesh with 100m grid resolution, thus doubling the grid points in each of the three coordinate directions. 
This re-gridding accomplished two objectives: 1) reduced the required computer resources (both CPU and disk 
space), and 2) accelerated the data post-processing. The re-gridding was based on the algorithms (Clark and Farley 
1984) that preserved accuracy and conservation of mass and energy. Comparison of interpolated fields between the 
200 m and 100m grids at the time of refinement verified that the interpolation retained the overall structure of the 
system. The re-gridded simulation was run for an additional 1 hr and 15mm to capture the evolution of the MCS in 
time. In both phases of the simulation, the model domain was allowed to translate with the MCS, thus relaxing the 
need for a much larger computational domain. 

The environmental conditions for TASS were initialized using a forecast sounding from the Weather Research 
and Forecast (WRF) model. The WRF outer most domain was bounded between 126.38°W and 69.72°W in the 
longitude and 21.59°N and 56.7 1°N in the latitude. The middle domain was bounded between 111.65°W and 
77.72°W in the longitude and 34.54°N and 54.07°N in the latitude. The innermost domain was bounded between 
100.02°W and 89.83°W in the longitude and 41.91°N and 48.30°N. The WRF outer-most domain had a horizontal 
grid resolution of 40 km in the horizontal and the two nested domains had horizontal resolutions of \3.33km and 
AAAkm respectively. Figure 5 shows the WRF computational domain with higher-resolution inner nests and the 
terrain heights for the inner-most nest. The simulation was initialized using North American Regional Reanalysis 
(NARR) data from the National Climate Data Center (NCDC) for 17 July 2008 at 1200UTC. See Mesinger et al. 
(2006) for details on the NARR model. 

The TASS input sounding (see Figure 6) was extracted from southern Minnesota for 1500UTC WRF forecast. 
The low levels of the sounding were modified to agree with observed cloud base heights. The environment may be 
characterized as unstable and moist with winds mostly from the west. The strongest winds were about 28 m/s (55 kts) 
and located near the tropopause. The freezing level was at 4100m AGL. No horizontal variation is assumed in the 
TASS initial conditions. Convection with cloud and precipitation development was triggered in TASS by imposing 
an artificial bubble of warm air. Three-dimensional fields evolved as convection interacted with the environmental 
fields. 
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Figure 4. Minneapolis NEXRAD (0.55° scan) image at 1523UTC. NEXRAD Station ID is KMPX. The white 
square shows the approximate location of the TASS 100km x 100km model domain. 
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Figure 5. The left panel shows the WRF computational domain with higher resolution nests and the right 
panel shows the terrain contours for the inner most computational domain. The terrain maximum value in 
the plot is 846.13 m and the minimum value is 175.97/w. 



-30 -20 -10 0 10 20 30 40 

Temperature (C) 


Figure 6. Skew-T chart of atmospheric sounding representing the environment for the Minnesota MCS. 
Shown are temperature, dewpoint, and horizontal wind velocity as a function of atmospheric pressure for 
southern Minnesota. 
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2. Validation and Analysis of Simulation 

A matching of the model time with UTC time shows that the simulation is best representative of conditions: 
corresponding UTC time = simulation time + 14:43; e.g., the TASS fine-grid results at 30 min corresponds to 
1514UTC. This was done with the aid of NEXRAD Level 2 data. A brief description of the NEXRAD products is 
given in Appendix A. In space, the TASS x-coordinate of x = 26,855 m corresponds to 93.77°W longitude, and the 
y-coordinate of y = 8,864 m corresponds to 44.22°N latitude. TASS simulation results were compared with the 
NEXRAD observations at Minneapolis (Table 2) and other surface observations. Similarities between the 
observations and simulation include a persistent MCS with many deep and active cells (updrafts) along its southern 
periphery. Other cells near center and northern regions can also be seen in both the observed data and simulation. 
Multiple updrafts fed the system with a large anvil spreading to the east. Hail developed in the simulation and was 
observed in the actual MCS at three different locations (see Table 3). The simulated and the observed storm tops 
obtained similar elevations of 15.5 km AGL (52,000 ft MSL). Maximum accumulated precipitation from the model 
was 39 .6mm (1.56m) and the maximum accumulated precipitation estimated from NEXRAD radar was between 
38 mm (1.5 in) and 63.5 mm (2.5 in). Observed Storm movement was around 17.5mA from the west and the simulated 
was about 16.5mA from the west. At storm mid-levels, a ridge of moderate reflectivity extends south eastward from 
the center of the system which can be seen on both the observed data and the simulation results. The simulated peak 
surface wind speed was approximately 30mA (67 mph) compared with the peak measured surface wind of 22.35 m/s 
(50 mph) recorded at Freeborn. The simulation differs from observations in system size. The observed system is 
larger than the modeled system (which was limited by the 100 km x 100 km domain). Also, the simulation did not 
generate cells along the eastern southern boundary that extend to the southeast along a cold front. A depiction of the 
simulated convective cloud is shown in Figures 7-8. Overshooting tops and broad eastward spreading anvil are 
apparent. Clusters of convective cells feeding the anvil are also apparent. 


Table 2: NEXRAD Station Data 


Station ID 

KMPX 

Location 

Minneapolis, MN 

Latitude 

44° 50’ 56” N 

Longitude 

93° 33’ 55” W 

Elevation (MSL) 

288 m (945/0 


Table 3: Model Comparison with Observations 


Variable 

TASS 

Observed 

Orientation of MCS Propagation 

West to East 

West to East 

Peak Storm Tops (AGL) 

15 .5km 

15.5 km 

Storm Motion (from) 

16.5 m/s (west) 

17.5mA (west) 

Maximum accumulated Precipitation 

39 .6mm (1.56z«) 

38 mm - 61 mm (1.5 in - 2 Ain)* 

Maximum Radar Reflectivity 

73dBZ 

6 8dBZ 

Peak Hail Size (diameter) 

14 mm 

9 mm continuous 

25 mm observed near Wasera between 
1430UTC-1500UTC 

18mm observed at Blue Earth (1536UTC) 
and Freeborn (1600UTC) 

Peak Surface Gust 

30 m/s 

22 m/s recorded at Freeborn (1600UTC) 
Trees reported down at other locations 


*Estimated from NEXRAD data. 


A composite of TASS RRF between the surface and 5 km is compared with the observed NEXRAD fields in 
Figures 9-10. The TASS radar reflectivity field (RRF) is computed within each grid cell from the size distribution 
of hydrometeors, but does not account for size and geometry of the radar beam. Both simulation and NEXRAD 
show the southwestern line of strong convective cells that is generated along the cooler downdrafts of air exiting the 
southern flank of the MCS. Cells near the center of the MCS are stronger and more persistent in TASS than those in 
the NEXRAD data. As noted earlier, the simulation does not capture the broken cells on the east side of the MCS. 
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3. Turbulence 


To define the turbulence characteristics, turbulence hazard fields were generated according to Eqs. (l)-(3) 
assuming an in-flight Boeing 737-300 with a weight of 100 klbs. The turbulence hazard is dependent upon aircraft 
type, with lighter aircraft generating stronger loads than heavy aircraft when encountering the same turbulence field. 
The simulated radar reflectivity factor and the turbulence intensity are shown at an elevation of 9Akm AGL 
(-32,000 ft MSL) in Figures 11-12. These fields are at a level slicing through the anvil, but well below the storm 
tops. First note, that radar reflectivity greater than 5 dBZ covers a broad area that is -50 km wide and almost 100 km 
long. Relatively small areas of radar reflectivity exceeding 30 dBZ are present, and are associated with vigorous 
convective cells. Turbulence with RMS normal load acceleration exceeding O.lg occurs within and in proximity to 
the higher reflectivity regions. However, spots with moderate to severe turbulence can be found within the low 
reflectivity regions of the MCS. Peak values of the RMS normal load acceleration at this altitude are in excess of 
0.4g. The figures also show the presence of many pockets of hazardous turbulence. Based on estimates appropriate 
for a Boeing 737-300 in-flight with a weight of 100 klbs, the normal accelerations due to turbulence would have 
exceeded 0.6g RMS of the normal load acceleration or 1.8g peak normal load in these pockets. 



Figure 7. TASS simulated cloud field at time = 15 min into the simulation. The view is from the south. Every 
other point in x- y- and z-direction is plotted. The simulation domain was bounded within a 100km x 100km x 
11km domain. Mesh resolution was set to Ax = Ay = Az = 100 m. 



Figure 8. TASS simulated cloud field at time = 15min into the simulation. The view is from the south. Every 
other point in x-y- and z-direction is plotted. The simulation domain was bounded within a 100km x 100km x 
11km domain. Mesh resolution was set to Ax = Ay = Az = 100 m. 
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Figure 9. Comparison of the KMPX NEXRAD 1.39° scan (left) with TASS low-level composite Radar 
Reflectivity Factor (right). Each panel depicts 100&/W x 100A:/w area. NEXRAD is at 1544UTC and TASS 
data is at time = 60 min into the simulation. 




Figure 10. Comparison of the KMPX NEXRAD 1.39° scan (left) with TASS low-level composite Radar 
Reflectivity Factor (right). Each panel depicts 100&/W x 100 km area. NEXRAD is at 1557UTC and TASS 
data is at time = ISmin into the simulation. 
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Figure 11. Radar reflectivity field and the root mean square of the normal load acceleration (a An ) computed 
from simulation at 60 min for an altitude of 9Akm AGL (~32,000/t MSL). The O.lg cr An contour is also 
depicted on radar reflectivity plot. Turbulence characteristics are based on data for an in-flight B737-300 
with a weight of 100 klbs. Axis coordinates are in meters. The simulation domain was bounded within a 
199km x 199km x 11km domain. Mesh resolution was set to Ax = Ay = Az = 100 m. 
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Figure 12. Radar reflectivity field and the root mean square of the normal load acceleration (a An ) computed 
from simulation at ISmin for an altitude of 9Akm AGL (~32,000/t MSL). The O.lg a An contour is also 
depicted on radar reflectivity plot. Turbulence characteristics are based on data for an in-flight B737-300 
with a weight of 100 klbs. Axis coordinates are in meters. The simulation domain was bounded within a 
199km x 199km x 11km domain. Mesh resolution was set to Ay = Ay = Az = 100 m. 
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B. Maritime Convection Induced Turbulence (3 August 2009) 

Isolated convective systems over ocean pose a serious turbulence hazard risk to aviation operations. Encounters 
with such systems are usually sudden and unexpected. In NASA’s Fall-2000 flight experiments, it was found that 
severe turbulence was often within overshooting tops and in growing cumulus turrets. Hamilton and Proctor (2002) 
reported that such encounters can occur in regions of low radar reflectivity values (peak values ranging from 8 dBZ 
to 33 dBZ). These small RRF values are difficult to detect by ground as well as airborne sensors. A likely example 
of such an encounter was reported on 3 August 2009, when the Continental Airlines Flight 128 encountered 
turbulence approximately 67 nautical miles north of Bavaro, Dominican Republic. The aircraft was a Boeing 767- 
224 flying at an altitude of 10. 91 km (36,000 ft). The NTSB incident report (NTSB ID: DCA09IA071) suggested that 
isolated and short-lived cumulus congestus to cumulonimbus clouds were developing near the flight path of the 
airplane. The flight had originated from Rio de Janeiro, Brazil and was bound for Houston, Texas but had to be 
diverted to Miami International Airport due to injuries to the passengers and crew members. Thirty three passengers 
and five flight attendants required medical attention as a result of this turbulence encounter. 

In this section, the simulation of a small and isolated tropical maritime convective event that developed south of 
Puerto Rico is examined. The storm formed off-shore and remained over water during its lifecycle. Throughout its 
lifecycle the storm remained within 60 nautical miles of the NEXRAD radar site at San Juan, Puerto Rico (station 
ID TJUA), and thus the event had reasonable coverage from ground-based radar. The NEXRAD station coordinates 
and elevation are given in Table 4. On 3 August 2009, scattered convective cells were observed by the San Juan 
NEXRAD radar in the vicinity of Puerto Rico. The streamlines at the 850 mb pressure level (Figure 13) from the 
North American Regional Reanalysis show that the flow was primarily easterly around 0900UTC. 

The convective-cell selected for TASS simulation was identified from the NEXRAD Level 3 data at 0819UTC. 
The storm started to weaken and dissipate around 0935UTC. The observed storm moved from the east and 
developed a slight tilt towards the north with time. The average speed of the storm was 6.88 m/s (13 Akts) at 96 
degrees. The observed echo tops and radar reflectivity at different elevation angles are shown in Figure 14. Echo 
tops as high as 12.3 km (>40,000//) were observed by the NEXRAD radar. No probability of hail was reported for 
the storm in the NEXRAD hail index. Radar reflectivity with values greater than 50 dBZ was observed at low levels, 
which is indicative of heavy rainfall at that time. A vertical cross section of the NEXRAD data is shown in Figure 
15. The details on storm trajectory, strength and precipitation are given in Table 5. All of the data in Table 5 was 
obtained from NEXRAD Level 3 output except the echo tops which, were estimated by analyzing the Level 2 
NEXRAD data. A description of the NEXRAD data in Table 5 is given in Appendix A. 



Figure 13. Streamlines at 850 mb pressure level on 3 August 2009 at 0900UTC. The data was obtained 
from the North American Regional Reanalysis (NARR) available from the National Climate Data Center 
(NCDC). 
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Figure 14. NEXRAD Level 2 data from the National Climate Data Center (NCDC). 3 August 2009 
(0916UTC). Echo tops (kft) are shown in the top left panel and the radar reflectivity (dBZ) at different 
elevation angles is shown in the other three panels. The island of Puerto Rico is located to the north of the 
cell. 



Figure 15. A vertical cross section of the radar reflectivity from the NEXRAD Level 2 data is shown in the 
figure. Source: National Climate Data Center (NCDC). 3 August 2009 (0916UTC). 
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Table 4: NEXRAD Station Data 


Station ID 

TJUA 

Location 

San Juan, PR 

Latitude 

18° 06’ 57” N 

Longitude 

66° 04’ 40” W 

Elevation (MSL) 

902 m (2958//) 


Table 5: Characteristics of the Storm - 3 August 2009 


Time 

(UTC) 

Lat 

(deg) 

Lon 

(deg) 

Range 

(nmi) 

Azimuth 

(deg) 

Echo Top 
(km) 

VIL 

(kg/m 2 ) 

Max 

Reflectivity 

(dBZ) 

One Hour 
Precip 
(mm) 

Storm 
Movement 
(deg | m/s) 

09 : 03 

17 . 659N 

66 . 810W 

50 

237 

8.6 

12 

50.5 

50 . 8 

93 

4.6 

09 : 10 

17 . 666N 

66 . 856W 

52 

239 

11.6 

13 

50.5 

50 . 8 

94 

5 . 1 

09 : 16 

17 . 671N 

66 . 879W 

53 

240 

12 .3 

12 

50.0 

50 . 8 

95 

5.6 

09:23 

17 . 682N 

66 . 926W 

55 

242 

12 .3 

10 

50 . 0 

50 . 8 

96 

6 . 1 

09:29 

17 . 681N 

66 . 965W 

57 

243 

11.6 

6 

48 . 0 

50 . 8 

96 

6.6 

09:35 

17 . 697N 

67 . 012W 

59 

245 

8 . 8 

3 

42 . 5 

38 . 1 

96 

6.6 


1. Simulation Setup 

The TASS computational domain was defined by a horizontal mesh of 18 km x 18 km, which extended vertically 
between the lower boundary representing the ocean’s surface to the model top set at 16.5 km. The grid resolution of 
the data set was 60 m in all coordinate directions (303 x 303 x 277 points). The model domain was allowed to 
translate with the simulated convective cell, thus relaxing the requirement for a larger computational domain. 

The environmental conditions for the TASS simulation were initialized using a forecast sounding from the WRF 
mesoscale weather prediction model. The WRF simulation domain consisted of an outer most domain bounded 
between 100.67°W and 31.43°W in the longitude and 5.21°S and 38.85°N with a mesh resolution of 66.5 km. Two 
higher resolution nests were defined within the outermost domain. Domain 2 was bounded between 85.9 1°W and 
46.2 1°W in the longitude and 3.8 1°N and 31.1°N in the latitude with a mesh resolution of 22.1 6km. The innermost 
domain was bounded between 76.37°W and 55.95°W in the longitude and 10.36°N and 25.05°N with a mesh 
resolution of 1.3%km. The WRF computational domain and the terrain for the innermost domain are shown in 
Figure 16. 

The WRF mesoscale simulation was initialized using Global Forecast System (GFS) data from the National 
Climate Data Center (NCDC) for 2 August 2009 at 1200 UTC. The simulation was run for 36 hours of forecast 
time. The output for the innermost domain was generated every 15 min. The sounding forecast (see Figure 17) was 
extracted from the WRF simulation high-resolution nest for 3 August at 0915 UTC. The low levels of the sounding 
were modified to agree with observed cloud base heights. The sounding indicates a moist environment favorable for 
single-cell convection with relatively small diameters. The initial environment assumed in TASS varies only in the 
vertical condition in accordance with the sounding in Figure 17. Convection with cloud and precipitation 
development was triggered in TASS by imposing the artificial bubble of warm air. Three-dimensional fields 
evolved as convection interacted with the environmental fields. Figure 18 shows the vertical cross sections of 
temperature and relative humidity through the location of extracted sounding (from the WRF simulation). A deep 
convection layer in the region surrounding the sounding location was conducive for triggering intense tropical 
thunderstorms. 
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Figure 16. The left panel shows the WRF computational domain with higher resolution nests and the right 
panel shows the terrain contours for the inner most computational domain. The terrain maximum value in 
the plot is 4447.95/w and the minimum value is 0 m. 
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Figure 17. Skew-T chart of atmospheric sounding representing the environment for the location southwest of 
Puerto Rico. Shown are temperature, dewpoint, and horizontal wind velocity as a function of atmospheric 
pressure. 
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Figure 18. WRF simulation data for 3 August 2009 at 0915UTC. Vertical cross-section of relative humidity 
and temperature through the location of the storm (latitude = 17.67°N and longitude = 66.86°W) are shown in 
the figure. 


2. Validation and Analysis of Simulation 

A matching of the model time with UTC time shows that the TASS simulation time = 28 min is best 
representative of the conditions at 0930UTC. This was done with the aid of NEXRAD Level 2 data. In space, the 
TASS x-coordinate of x = -11,151m is equivalent to 66.93°W longitude and, the y-coordinate of y = 1024m is 
equivalent to 17.67°N latitude. The TASS simulation results are compared with observations in Table 6. The 
simulated storm motion (both speed and direction) closely match the observed system. The simulated echo tops are 
also in good agreement with the observations. The computed radar reflectivity however is larger than observed from 
NEXRAD. In Figures 19-20, the TASS simulated cloud fields are shown at two different times (20 min and 32 min). 
The interaction of the convective system with the strong environmental shear at higher altitudes can be seen in 
Figure 20. Figure 21 shows a cross section of the TASS computed radar reflectivity in x-z plane for y = 2071m at 
time = 30 min. 


Table 6: Model Comparison with Observations 


Variable 

TASS 

Observed 

Orientation of Storm 

East - West (slight tilt towards North) 

East - West (slight tilt towards North) 

Peak Storm Tops (AGL) 

12.5 km 

12.3 km 

Storm Motion 

6.1 m/s 95 degrees 

6.8 m/s 96 degrees (average values) 

Max Radar Reflectivity 

6(klBZ 

50.5dBZ 


18 

American Institute of Aeronautics and Astronautics 



Figure 19. TASS simulated cloud field at time = 20 min into the simulation. View is from the south- 
southwest. Every other point in x- y- and ^-direction is plotted. The simulation domain was bounded within 
an 18km x 18km x 16.5 km domain. Mesh resolution was set to Ax = Ay = Az = 60 m. 



Figure 20. TASS simulated cloud field at time = 32 min into the simulation. View is from the south- 
southwest. Every other point in x- y- and z-direction is plotted. The simulation domain was bounded within 
an 18 km x \8km x 16.5 km domain. Mesh resolution was set to Ax = Ay = Az = 60 m. 


19 

American Institute of Aeronautics and Astronautics 




15000 


10000 

g 

N 


5000 



-20000 


-15000 -100DG 

x{m) 


-5000 


Figure 21. Vertical cross section (. x-z plane at y = 2611m) of the TASS simulated radar reflectivity field at 
time = 30 min into the simulation. Every other point in x-y- and ^-direction is plotted. The simulation domain 
was bounded within an 18km x 18 km x 16.5km domain. Mesh resolution was set to Ax = Ay = Az = 60 m. 


3. Turbulence 

The turbulence characteristics for this case are shown in Figure 22 along with the radar reflectivity factor. The 
root mean square of the normal load acceleration is computed assuming an in-flight Boeing 767 with a weight of 
300 klbs. Light to moderate turbulence intensities occur within the top levels of this storm and coincide with low 
values of radar reflectivity. Convection induced turbulence was encountered by a Boeing 767 at an elevation of 
\0.91km (36,000 ft) in a similar environment that was northeast of Puerto Rico on this date (3 August 2009). 
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Figure 22. Radar reflectivity field and the root mean square of the normal load acceleration (GAn) 
computed from simulation at time = 30 min for an altitude of 16.98km AGL. Turbulence characteristics are 
based on data for an in-flight B-767 with a weight of 300 klbs. Every other point in x- y- and z-direction is 
plotted. The simulation domain was bounded within an 18km x 18km x 16.5km domain. Mesh resolution was 
set to Ay = Ay = Az = 60/;/. 
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C. Indiana Hailstorm (19 June 2009) 

A supercell convective storm formed on the border of Fulton County and Marshall County on 19 June 2009 at 
approximately 2000UTC. The storm was located southwest of the NEXRAD radar site at Fort Wayne, Indiana. The 
prevalent synoptic conditions at 2100UTC on 19 June 2009 are shown in Figure 23 along with the infrared imagery 
from the GOES-E satellite. The storm persisted for a couple of hours and moved westnorthwest (278 degrees) at a 
speed of 10.7mA (21 kts). Throughout its lifecycle the storm remained within 60 km of the NEXRAD radar site at 
Fort Wayne, Indiana (station ID KIWX). The NEXRAD station coordinates and elevation are given in Table 7. The 
storm characteristics (storm location, trajectory and strength) obtained from an analysis of the NEXRAD L2 and L3 
data are summarized in Table 8. 

The supercell resulted in a hailstorm that produced softball sized hail (114.3mm) and generated a wall cloud. 
Observations of hail were reported from 2000UTC until 2206UTC. Hail ranging from pea to grapefruit size and 
lasting for 20mm was reported in Fulton County at 2030UTC. Possibility of hail (49.2mm) was reported by the 
NEXRAD data in the Marshall County at 2016UTC. Hail (possible size ranging from 54.6mm to 60.9mm) was also 
reported in Kosciusko County at 2055UTC and again at 2118UTC, as the storm moved east. No substantial wind 
damage was reported for Fulton County and the storm in general had low shear in the lower 1 km of the boundary 
layer. 

Analysis of the NEXRAD data indicated radar reflectivity in excess of lOdBZ and storm tops of \lkm 
(>55,000 ft). The echo tops and the radar reflectivity values (from Level 2 data) at different elevation angles are 
shown in Figure 24 for 2044UTC. Radar reflectivity of 65 dBZ at low levels was observed, indicating very heavy 
rainfall and the possibility of marble to golf ball sized hail at that time. The maximum echo tops recorded at this 
time were higher than ll.lkm (56,000 ft). A vertical cross section of the Level 2 NEXRAD data is shown in Figure 
25. A region of reflectivity with values as high as lOdBZ can be seen in Figure 25. A detailed description of the 
storm is given by B. J. Simpson of the National Weather Service (NWS) in a presentation posted on the NWS 
website (Simpson 2009). Simpson also discusses the reasons why this particular supercell with sufficient energy 
(MLCAPE in excess of 4000 J/kg) did not produce tornados. One of the reasons for the failure of tomadogenesis is 
attributed to weak low level shear. 



Figure 23. Surface weather map for the continental United States overlaid with GOES-E IR imagery at 
2100UTC on 19 June 2009. Source: National Centers for Environmental Prediction (NCEP). 
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Figure 24. NEXRAD Level 2 data from the National Climate Data Center (NCDC). 19 June 2009 
(2044UTC). Echo tops ( kft ) are shown in the top left panel and the radar reflectivity ( dBZ) at different 
elevation angles is shown in the other three panels. The white diamond represents the observed hail size of 
87.4/w/w. 



Figure 25. NEXRAD Level 2 data from the National Climate Data Center (NCDC). 19 June 2009 
(2044UTC). A vertical cross section of the radar reflectivity is shown in the figure. 
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Table 7: NEXRAD Station Data 


Station ID 

KIWX 

Location 

Fort Wayne, IN 

Latitude 

41° 21’ 32” N 

Longitude 

85° 42’ 00” W 

Elevation (MSL) 

321.6m (1055/0 


Table 8: Characteristics of the Storm - 19 June 2009 


Time 

(UTC) 

Lat 

(deg) 

Lon 

(deg) 

Range 

(nmi) 

Azimuth 

(deg) 

Echo Top 
(km) 

VIL 

(kg/m 2 ) 

Max 

Reflectivity 

(dBZ) 

One 

Hour 

Precip 

(mm) 

Storm 
Movement 
(deg | m/s) 

20:01 

41 . 172N 

86 .410W 

34 

251 

13.5 

36 

59 . 0 

2 . 54 

254 | 

11 . 7 

20:06 

41 . 178N 

86 .399W 

33 

251 

14 . 5 

50 

59 . 0 

6.35 

256 | 

11.2 

20 : 11 

41 . 172N 

86 .339W 

31 

249 

15.2 

62 

65 . 0 

12 . 7 

260 | 

11.2 

20:16 

41 . 176N 

86 .294W 

29 

248 

15 . 0 

65 

63 . 0 

12 . 7 

263 | 

11.2 

20:21 

41 . 167N 

86 .240W 

27 

245 

15.4 

71 

64 . 0 

19 . 1 

265 | 

11 . 7 

20:25 

41 . 154N 

86 .207W 

26 

242 

15.3 

68 

66 . 0 

19 . 1 

269 | 

11 . 7 

20:30 

41 . 128N 

86 . 186W 

26 

238 

15 . 9 

81 

69 . 0 

25.4 

273 | 

11 . 7 

20:34 

41 . 119N 

86 . 152W 

25 

235 

17.2 

80 

66 

25.4 

280 | 

11.2 

20:39 

41 . 117N 

86 . 123W 

24 

233 

16 . 8 

78 

67 

25.4 

284 | 

11.2 

20:44 

41 . 102N 

86 . 077W 

23 

228 

17.2 

78 

66 

25.4 

287 | 

11 . 7 

20:48 

41 . 113N 

86 . 061W 

22 

228 

16 . 1 

74 

68 

25.4 

289 | 

11.2 

20 : 53 

41 . 111N 

86 . 027W 

21 

225 

15.3 

72 

65 . 0 

31 . 8 

289 | 

10 . 7 

20 : 57 

41 . 107N 

85.989W 

20 

221 

15.3 

68 

66 . 0 

38 . 1 

289 | 

10.2 

21 : 02 

41 . 079N 

85 . 979W 

21 

217 

14 .4 

62 

64 . 0 

38 . 1 

288 

9.7 


1. Simulation Setup 

The TASS computational domain was bounded in the horizontal within a mesh of 10km x 10km. The vertical 
extant was set to 18.5 km. The grid resolution of the data set was 100m in all coordinate directions (703 x 703 x 188 
points). The TASS simulation was initialized using a sounding from the Weather Research and Forecast model 
simulation. The WRF simulation domain consisted of an outer most domain bounded between 122.92°W and 
46.07°W in the longitude and 20.79°N and 58.14°N with a mesh resolution of 18.0 km. Two higher resolution nests 
were defined within the outermost domain. Domain 2 was bounded between 107.82°W and 65.03°W in the 
longitude and 28.82°N and 51.62°N in the latitude with a mesh resolution of 6.0km. The innermost domain was 
bounded between 95.95°W and 76.65°W in the longitude and 35.31°N and 46.19°N with a mesh resolution of 
2.0km. High-resolution (1 km) terrain and land use datasets were used for the innermost domain and coarse 
resolution datasets were used for outer domains. The WRF computational domain and the terrain for the innermost 
domain are shown in Figure 26. The simulation was initialized using North American Regional Reanalysis (NARR) 


23 

American Institute of Aeronautics and Astronautics 


data from the National Climate Data Center (NCDC) for 19 June 2009 at 1200 UTC. The simulation was run for 24 
hours. The output for the innermost domain was generated every 30 min. 

The sounding forecast (see Figure 27) was extracted from the WRF simulation high-resolution nest for 19 June 
at 2000UTC. The low levels of the sounding were modified to agree with observed cloud base heights. The initial 
environment assumed in TASS varies only in the vertical condition in accordance to the sounding in Figure 27. 
Convection with cloud and precipitation development was triggered in TASS by imposing an artificial bubble of 
warm air. Three-dimensional velocity and thermodynamic quantities evolved as convection interacted with the 
environmental fields. 



Figure 26. The left panel shows the WRF computational domain with higher resolution nests and the right 
panel shows the terrain contours for the inner most computational domain. The terrain maximum value in 
the plot is 1726.62/w and the minimum value is 15.47/w. 
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Figure 27. Skew-T chart of atmospheric sounding representing the ambient environment. Shown are 
temperature, dewpoint, and horizontal wind velocity as a function of atmospheric pressure. 
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2. Validation and Analysis of Simulation 

A matching of the model time with UTC time shows that the TASS simulation time = 53 min is best 
representative of the conditions at 2030UTC. In space, the TASS x-coordinate of x = 25 km corresponds to 
86.178°W longitude and, the y-coordinate of y= \2.15km corresponds to 41.09°N latitude. The temporal and spatial 
mapping between the TASS and the observed data was done with the help of NEXRAD radar data. The TASS 
simulation results are compared with observations in Table 9. In Figure 28, the TASS simulated cloud field is 
shown at time = 50 min and 15s into the simulation. The TASS computed RRF fields at two different heights are 
shown in Figure 29. The formation of the hook echo and the gust front can be seen in the figure. These features 
were also observed in the actual hailstorm. 


Table 9: Model Comparison with Observations 


Variable 

TASS 

Observed 

Orientation of Storm 

West - East 

West - East (slight tilt to the south) 

Peak Storm Tops (AGL) 

1 6km 

\12km 

Storm Motion 

ll.lm/s 255 degrees 

\0.1m/s 278 degrees (average values) 

Maximum Radar Reflectivity 

IHAdBZ 

6 9.0dBZ 



Figure 28. TASS simulated cloud field at time = 50min:\5s. View is from the southwest. The simulation 
domain is bounded within a 10km x 10km x 18.5km domain. Mesh resolution was set to Ax = Ay = Az = 100 m. 
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Figure 29. The TASS computed RRF field in the x-y plane at an altitude of 2200 m is shown in the left panel 
and at an altitude of 6000 m is shown in the right panel. Time = 53 min into the simulation. 
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3. Turbulence 

The RMS of the normal load acceleration (cr An ) was computed from the TASS simulation data assuming an in- 
flight Boeing B737-300 with a weight of 100 klbs. Significant cr An values (cr An > 0.2g) occur within the top levels of 
this storm and coincide with low values of radar reflectivity. These values correspond to peak normal loads (]A«|) 
greater than 0.5g. The 0.2 g isosurface of the RMS normal load acceleration from TASS simulation is shown in 
Figure 30 along with the 35 dBZ radar reflectivity isosurface. The maximum cr An value at this time was 0.77 g. The 
TASS computed RRF and the cr An fields in the x-y plane at an altitude of 10 km are shown in Figure 31. 
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Figure 30. Radar Reflectivity Factor (RRF) isosurface (35 dBZ) is shown in green and the 0.2 g a An 
isosurface is shown in red. A Boeing 737-300 aircraft with a weight of 100 klbs is assumed in a An calculation. 
Every point in ^-direction is plotted and every other point in x- and j-direction is plotted. The simulation 
domain was bounded within a 10km x 10km x 18.5 km domain. Mesh resolution was set to Ax = Ay = Az = 
100 m. The TASS simulation time = 53 min. 
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Figure 31. The TASS computed Radar Reflectivity Factor field in the x-y plane at an altitude of 10000/w is 
shown in the left panel and the a An field on the same plane is shown in the right panel. Every other point in x- 
and j-direction is plotted. A Boeing 737-300 aircraft with a weight of lOOklbs is assumed in the a An 
calculation. Time = 53min into the simulation. 
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D. Puerto Rico (17 July 2010) 

On 17 July 2010, a strong gust front formed on the eastern shore of Puerto Rico and started moving 
westnorthwest. The 1 km resolution visible imagery (Figure 32) from GOES- 13 shows the outflow boundary as the 
front propagated westwards. The Intertropical Convergence Zone (ITCZ) is normally south of Puerto Rico in July 
and it is possible that the initial forcing for the gust front may have been provided by convection due to the 
convergence in the ITCZ. There were unconfirmed reports of a waterspout/tomado near Mayaguez around 
1645UTC along with reports of damage to trees and structures. Multiple injuries were also reported and surface 
wind gusts in excess of 22.35 m/s (50 mph) were recorded at various locations in the path of the gust front. An 

enhanced-V signature (AT 10°C) can be seen in the MODIS infrared imagery (Figure 32). The enhanced-V 

signature is formed due to winds going around overshooting tops and is indicative of severe thunderstorms (McCann 
1983). 

The thunderstorm event over northwestern Puerto Rico was detected at 1630UTC by the NEXRAD radar site 
near San Juan, Puerto Rico (NEXRAD station ID TJUA). Several convective storms with multiple cells formed in 
the region. These cells split and merged into each other as the gust front moved over northwest Puerto Rico. The 
storm moved briskly at approximately 1 1.82mA towards the westnorthwest. The storm was later observed to move 
out over the ocean and dissipate by 1800UTC. Hail (size ranging from 21. 6mm to 24.1mm) was estimated from the 
NEXRAD at 1706UTC. Maximum radar reflectivity of 51dBZ and echo tops as high as \43km (4!kft) were 
observed during the storm lifetime. The radar reflectivity (Level 2 data) and echo tops at 1706UTC at different 
elevation angles is shown in Figure 33. A vertical cross section of the Level 2 NEXRAD data is shown in Figure 
34. The details on storm trajectory, strength and precipitation are given in Table 10. 




Figure 32. Visible imagery from NOAA’s GOES-13 satellite for 17 July 2010. 1655UTC (top left), 

1745UTC (top right), 1845UTC (bottom left). The Aqua MODIS IR imagery at 1744UTC is shown in the 
bottom right panel. Source: Cooperative Institute for Meteorological Satellite Studies (CIMSS). 
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Figure 33. NEXRAD Level 2 data from the National Climate Data Center (NCDC). 17 July 2010 
(1706UTC). Echo tops (kft) are shown in the top left panel and the radar reflectivity (dBZ) at different 
elevation angles is shown in the other three panels. 



Figure 34. NEXRAD Level 2 data from the National Climate Data Center (NCDC). 17 July 2010 
(1706UTC). A vertical cross section of the radar reflectivity is shown in the figure. 
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Table 10: Characteristics of the Storm - 17 July 2010 


Time 

(UTC) 

Lat 

(deg) 

Lon 

(deg) 

Range 

(nmi) 

Azimuth 

(deg) 

Echo Top 
(km) 

VIL 

(kg/m 2 ) 

Max 

Reflectivity 

(dBZ) 

One Hour 
Precip 
(mm) 

Storm 
Movement 
(deg | m/s) 

16:31 

18 .344N 

66 . 866W 

47 

287 

9.2 

14 

51 

n/ a 

New 

16:42 

18 .392N 

66 . 923W 

51 

289 

12 . 1 

7 

52 

n/ a 

123 

1 

13 .3 

16:48 

18 .412N 

66 . 935W 

52 

290 

13 . 5 

5 

55 

n/ a 

127 

1 

11 . 7 

16 : 54 

18 .429N 

66 . 984W 

55 

290 

14 .4 

28 

55 

n/ a 

124 

1 

11 . 7 

17 : 00 

18 .429N 

66 . 984W 

55 

290 

14 . 0 

31 

53 

n/ a 

125 

1 

10.2 

17 : 06 

18 .440N 

67 . 017W 

57 

290 

14 . 1 

33 

56 

n/ a 

123 

1 

10.2 

17 : 11 

18 .467N 

67 . 044W 

59 

291 

14 . 1 

27 

56 

n/ a 

123 

1 

13.3 

17 : 24 

18 . 525N 

67 . 095W 

63 

293 

14 . 0 

26 

52 

76.2 

126 

1 

9 . 7 

17:30 

18 . 544N 

67 . 144W 

66 

293 

13 . 0 

26 

50 

76.2 

127 

1 

10.2 

17:35 

18 . 582N 

67 . 184W 

69 

294 

12 . 1 

26 

51 

50 . 8 

123 

1 

13 .3 

17:41 

18.596N 

67.216W 

71 

294 

13 . 1 

21 

53 

50 . 8 

129 

1 

11.2 

17:47 

18 . 602N 

67.232W 

72 

294 

13 .7 

18 

52.5 

50 . 8 

130 

1 

11.2 

17 : 53 

18 . 629N 

67.239W 

73 

295 

11.2 

75 

51 

50 . 8 

131 

1 

11.2 

17 : 59 

18 . 733N 

67.303W 

79 

298 

11.5 

20 

54 

50 . 8 

135 

1 

12.2 


1. Simulation Setup 

The TASS computational domain was bounded in the horizontal within a mesh of 45 km x 45 km. The model top 
was set to \1 2km. The grid resolution was set to 80m in all coordinate directions (566 x 566 x 216 points). The 
model domain was allowed to translate with simulated convective cell, thus relaxing the requirement for a larger 
computational domain. The TASS simulation was initialized from a sounding obtained from a WRF simulation. 

The WRF simulation domain consisted of an outer most domain bounded between 108.98°W and 23.13°W in the 
longitude and 9.86°S and 43.04°N with a mesh resolution of 3 6km. Two higher resolution nests were defined within 
the outermost domain. Domain 2 was bounded between 86.54°W and 45.57°W in the longitude and 3.90°N and 
31.27°N in the latitude with a mesh resolution of 12 km. The innermost domain was bounded between 75.56°W and 
56.54°W in the longitude and 11.52°N and 24.39°N with a mesh resolution of 4 km. High-resolution (1 km) terrain 
and land use datasets were used for the innermost domain. The WRF computational domain and the terrain for the 
innermost domain are shown in Figure 35. The simulation was initialized using Global Forecast System (GFS) data 
from the National Climate Data Center (NCDC) for 17 July 2010 at 1200UTC. The simulation was run for 12 
hours. The output for the innermost domain was generated every 15 min. The sounding forecast (see Figure 36) was 
extracted from the WRF high-resolution nest for 17 July at 1530UTC. The initial environment assumed in TASS 
varies only in the vertical condition in accordance to the sounding in Figure 36. Convection with cloud and 
precipitation development was triggered in TASS by imposing an artificial bubble of warm air. Three-dimensional 
velocity and thermodynamic quantities evolved as convection interacted with the environmental fields. 
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Figure 35. The left panel shows the WRF computational domain with higher resolution nests and the right 
panel shows the terrain contours for the inner most computational domain. The terrain maximum value in 
the plot is 2400.94 m and the minimum value is 0 m. 
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Figure 36. Skew-T chart of atmospheric sounding representing the ambient environment at the time of TASS 
initialization. Shown are temperature, dewpoint, and horizontal wind velocity as a function of atmospheric 
pressure. 
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2. Validation and Analysis of Simulation 

A matching of the model time with the UTC time shows that the TASS simulation at time = 40 min is best 
representative of conditions at 1712UTC. This was done with the aid of NEXRAD Level 2 data. The TASS 
coordinates of x = -23201m and y = 10698 m are equivalent to 67.065°W longitude and 18.482°N latitude. The 
TASS simulation captured one of the 3 or 4 storms that stretched from the southeast to the northwest. Simulation 
results are compared with the observations in Table 11. The model is able to capture the salient features of the 
storm such as peak echo tops and the storm motion and direction, as well as the shearing of the cloud system from 
the northwest to southeast. The angle of the shearing is more southward in the NEXRAD data. Also, the peak radar 
reflectivity is higher in the computations compared to the NEXRAD data. 

The TASS simulated cloud field at time = 43 min and 15s is shown in Figure 37. The tilt of the cloud at higher 
altitudes due to shear in the environment is apparent in the figure. The TASS computed radar reflectivity is 
compared with the San Juan NEXRAD 0.47° scan in Figure 38. 


Table 11: Model Comparison with Observations 


Variable 

TASS 

Observed 

Orientation of Storm 

Southeast - Northwest 

Southeast - Northwest 

Peak Storm Tops (AGL) 

13.72km 

14.38km 

Storm Motion 

I2.4m/s 115.14 degrees 

11.5 m/s 127 degrees (average values) 

Maximum Radar Reflectivity 

64.26dBZ 

56.(klBZ 



Figure 37. TASS simulated cloud field at time = 43min and 15s. View is from the south. This is around the 
time of maximum storm intensity. Every other point in x-, y- and z-direction is plotted. The simulation 
domain was bounded within a 45km x 45km x 11.2km domain. Mesh resolution was set to Ax = Aj = Az = 80 m. 
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Figure 38. Comparison of the TJUA NEXRAD 0.47° scan (left) with TASS RRF (right) at an altitude of 
1600 m. The white box has an area of 45km x 45km which is the size of the TASS computational domain. 
NEXRAD is at 1712UTC and TASS data is at 40 min into the simulation. 


3. Turbulence 

The RMS of the normal load acceleration (cr An ) was computed from the TASS simulation data assuming an in- 
flight Boeing B737-300 with a weight of 100 klbs. Significant cr An values (cr^ > 0.2g) occur within the top levels of 
this storm and coincide with low values of radar reflectivity. These values correspond to peak normal loads greater 
than 0.5g (see Figure 1). The maximum cr An value at this time in simulation was 0.46g. The 0.2 g <r An isosurface 
from TASS simulation is shown in Figure 39 along with the 35 dBZ radar reflectivity isosurface. A portion of the 
vertical velocity slice in the x-z plane coinciding with the 0.2g cr An isosurface is also shown in the figure. 
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Figure 39. Radar Reflectivity Factor (RRF) isosurface (35 dBZ) is shown in green and the 0.2g a An 
isosurface is shown in red. A portion of the vertical velocity (m/s) slice through the x-z plane at y = 12260 m is 
also shown in the figure. A Boeing 737-300 aircraft with a weight of 100 klbs is assumed in a An calculation. 
The TASS simulation time = 43 min and 15s. 
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E. Lincoln, Illinois (10 August 2010) 

An isolated storm was identified by the Lincoln, Illinois NEXRAD station (station ID KILX) at 21 17UTC on 10 
August 2010. The NEXRAD station data is given in Table 12. The storm was assigned the identifier B2. The 
storm was re-assigned the identifier OO at 2150UTC. The storm OO dissipated by 2245UTC. The storm remained 
north of a front that stretched through the middle of Illinois (Figure 40). It meandered in the region at a slow speed 
(~6kts) towards the east and remained in the Ford County. Throughout its evolution it remained within a 60 nautical 
miles radius around the radar site at Lincoln. The storm reached its maximum intensity around 2217UTC. The 
storm split into two cells and started to weaken and dissipate around 2245UTC. The cell completely dissipated by 
2308UTC. 

Probability of hail (possible sizes ranging from 19.1mm to 38.1mm) was reported for both the storms B2 and OO 
in the NEXRAD Level 2 data and the Level 3 hail index. The radar reflectivity (Level 2 data) at 2150UTC at 
different elevation angles are shown in Figure 41 . Radar reflectivity of 52.5 dBZ at low levels was observed which is 
indicative of heavy rainfall at that time. The maximum echo top at this time was \2.6km (41.3 kft). A vertical cross 
section of the Level 2 NEXRAD data is shown in Figure 42. A region of reflectivity with values in the range of 
5 6dBZ can be seen in Figure 42. 

The details on storm trajectory, strength and precipitation are given in Table 13. All of the data in Table 13 was 
obtained from NEXRAD Level 3 output except the echo tops which, were estimated by analyzing the Level 2 
NEXRAD data. The actual event, as detected by NEXRAD radar, lasted for about 2 hours, had storm tops upwards 
of 12.9 km (42.3 kft), and moved relatively slowly (~3 m/s) towards the east. 



Figure 40. Surface weather map for the continental United States overlaid with GOES-E IR imagery at 
2115UTC on 10 August 2010. Source: National Centers for Environmental Prediction (NCEP). 
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Figure 41. NEXRAD Level 2 data from the National Climate Data Center (NCDC). 10 August 2010 
(2150UTC). Echo tops (lift) are shown in the top left panel and the radar reflectivity ( dBZ) at different 
elevation angles is shown in the other three panels. The storm of interest is near Paxton, Illinois. 



Figure 42. A vertical cross section of the NEXRAD Level 2 data radar reflectivity is shown. 10 August 
2010 (2150UTC). Source: National Climate Data Center (NCDC). 
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Table 12: NEXRAD Station Data 


Station ID 

KILX 

Location 

Lincoln, IL 

Latitude 

40° 08’ 59” N 

Longitude 

89° 20’ 13” W 

Elevation (MSL) 

222.5m (730 ft) 


Table 13: Characteristics of the Storm - 10 August 2010 


Time 

(UTC) 

Lat 

(deg) 

Lon 

(deg) 

Range 

(nmi) 

Azimuth 

(deg) 

Echo 

Top 

(km) 

VIL 

(kg/m 2 ) 

Max 

Reflectivity 

(dBZ) 

One Hour 
Precip 
(mm) 

Storm 
Movement 
(deg | m/s) 

B2/21 : 17 

40 . 526N 

88 .315W 

52 

64 

8.1 

3 

52 . 0 

2.5 

new 

B2/21 :22 

40 . 512N 

88.307W 

52 

65 

8.1 

6 

54 . 0 

2.5 

312 

1 7.1 

B2/21 : 27 

40 . 505N 

88 .279W 

53 

66 

8.4 

9 

52.5 

6.3 

296 

5.6 

B2/21 : 36 

40 . 512N 

88 .307W 

52 

65 

n/ a 

16 

50 . 5 

6.3 

5 1 

1 . 0 

00/21:50 

40 . 524N 

88 .363W 

50 

63 

12.6 

22 

50 . 5 

12 . 7 

New 

00/21:55 

40 . 524N 

88 .363W 

50 

63 

12.2 

34 

53 . 0 

12 . 7 

163 

1.5 

00/22:00 

40 . 537N 

88 .372W 

50 

62 

12.6 

46 

58 . 5 

12 . 7 

160 

2.6 

00/22:06 

40 . 545N 

88 .352W 

51 

62 

12 . 8 

46 

58 . 5 

25.4 

181 

2.6 

00 / 22 : 12 

40 . 532N 

88 .344W 

51 

63 

12 . 9 

31 

59 . 0 

31 . 8 

208 

2.6 

00 / 22 : 17 

40 . 532N 

88 .344W 

51 

63 

12.5 

34 

57.5 

38 . 1 

206 

1 1-5 

00 / 22:23 

40 . 532N 

88 .344W 

51 

63 

12 .4 

27 

56.5 

38 . 1 

n/ a 

00 / 22:34 

40 . 526N 

88 .315W 

52 

64 

11.6 

12 

53 . 5 

n/ a 

n/ a 

00 / 22:40 

40 . 519N 

88 .287W 

53 

65 

9 . 1 

16 

51 . 0 

44 . 5 

270 

1 1-5 

00 / 22:45 

40 . 512N 

88 .307W 

52 

65 

8.9 

11 

54 . 0 

44 . 5 

285 

| 2.0 


1. Simulation Setup 

The TASS computational domain was bounded in the horizontal within a mesh of 42 km x Allan. The vertical 
extant was set to 18 km. The grid resolution was set to 80m in all coordinate directions (528 x 528 x 229 points). 
The model domain translates with simulated convective cell, thus relaxing the requirement for a larger 
computational domain. The TASS simulation was initialized from the WRF sounding. 

The WRF simulation domain consisted of an outer most domain bounded between 130.01°W and 50.8 1°W in the 
longitude and 20.33°N and 58.44°N with a mesh resolution of 60.0 km. Two higher resolution nests were defined 
within the outermost domain. Domain 2 was bounded between 109.12°W and 71.69°W in the longitude and 
30.63°N and 50.25°N in the latitude with a mesh resolution of 20.0 km. The innermost domain was bounded 
between 99.01°W and 81.55°W in the longitude and 36.28°N and 45.31°N with a mesh resolution of 6.66km. High- 
resolution (1 km) terrain and land use datasets were used for the innermost domain. The WRF computational domain 
and the terrain for the innermost domain are shown in Figure 43. The simulation was initialized using Global 
Forecast System (GFS) data from the National Climate Data Center (NCDC) for 10 August 2010 at 1200UTC. The 
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simulation was run for 24 hours. The output for the innermost domain was generated every 15 min. The sounding 
forecast (see Figure 44) was extracted from the WRF simulation high-resolution nest for 10 August at 2130UTC. 
The low levels of the sounding were modified to agree with observed cloud base heights. The initial conditions in 
TASS vary only in the vertical direction in accordance to the sounding in Figure 44. Convection with cloud and 
precipitation development was triggered in TASS by imposing an artificial bubble of warm air. Three-dimensional 
fields evolved as convection interacted with the environmental fields. 



110°W 100°W 90°W 80°W 70°W 0 100 200 300 400 500 600 700 800 900 1000 


Figure 43. The left panel shows the WRF computational domain with higher resolution nests and the right 
panel shows the terrain contours for the inner most computational domain. The terrain maximum value in 
the plot is 869.96 m and the minimum value is 82.98 m. 
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Figure 44. Skew-T chart of atmospheric sounding representing the ambient environment used for TASS 
initialization. Shown are temperature, dewpoint, and horizontal wind velocity as a function of atmospheric 
pressure. 
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2. Validation and Analysis of Simulation 

A matching of the model time with the UTC time shows that the TASS simulation at time = 37 min is best 
representative of conditions at 2217UTC. This was done with the aid of NEXRAD Level 2 data. The TASS 
coordinates of x = 8521.83m and y = 501.51m are equivalent to 88.34°W longitude and 40.536°N latitude. The 
TASS simulation results are compared with observations in Table 14. The TASS computed radar reflectivity is 
over-predicted but other measures such as, storm orientation, storm motion and peak echo tops are in good 
agreement with observations. 

In Figure 45, the TASS simulated cloud field is shown at time = 37mm into the simulation. In Figures 46-47, the 
cross sections of the TASS computed radar reflectivity in x-y plane at two different heights are compared with 
NEXRAD data at time = 37mm into the simulation. 


Table 14: Model Comparison with Observations 


Variable 

TASS 

Observed 

Orientation of Storm 

East - West 

East - West 

Peak Storm Tops (AGL) 

14.5 km 

12.9 km 

Storm Motion 

3.8 m/s 

2.8 m/s 

Maximum Radar Reflectivity 

16MBZ 

59.(klBZ 



Figure 45. TASS simulated cloud field at time = 37 min is shown. The view is from the south. Every point 
in z-direction is plotted and every other point in x- and j-direction is plotted. The simulation domain was 
bounded within a 42 km x 42 km x 18km domain. Mesh resolution was set to Av = Ay = Az = 80 m. 
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Figure 46. Comparison of the KILX NEXRAD 3.37° scan (left) with TASS RRF (right). The white box is 
42 km x 42 km in area and represents the size of TASS simulation domain. TASS x-y plane at an altitude of 
6720/w is shown in the figure. NEXRAD is at 2217UTC and TASS data is at 37 min into the simulation. Every 
other point in x- and y-direction is plotted. The simulation domain was bounded within a 42km x 42km x 
18km domain. Mesh resolution is set to Ax = Ay = Az = 80 m. 
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Figure 47. Comparison of the KILX NEXRAD 5.99° scan (left) with TASS RRF (right). The white box is 
42 km x 42 km in area and represents the size of TASS simulation domain. TASS x-y plane at an altitude of 
12000 m is shown in the figure. NEXRAD is at 2217UTC and TASS data is at 31min into the simulation. 
Every other point in x- and y-direction is plotted. The simulation domain was bounded within a 42 km x 42km 
x 18 km domain. Mesh resolution is set to Av = Ay = Az = 80 m. 
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3. Turbulence 

The RMS of the normal load acceleration (cr An ) was computed from the TASS simulation data assuming an in- 
flight Boeing B737-300 with a weight of 100 klbs. Significant cr An values (cr An > 0.2g) occur within the top levels of 
this storm and coincide with low values of radar reflectivity. These values correspond to peak normal loads greater 
than 0.5g. The 0.2g isosurface of the RMS normal load acceleration from TASS simulation is shown in Figure 48 
along with 35 dBZ radar reflectivity isosurface. The maximum cr An value at this time was 0.6g. The TASS computed 
RRF and the cr An fields in the x-y plane at an altitude of 1 1 120m are shown in Figure 49. 
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Figure 48. TASS computed Radar Reflectivity Factor isosurface of 35 dBZ is shown in green and the 0.2g 
a An isosurface is shown in red. A Boeing 737-300 aircraft with a weight of \Wklbs is assumed in a An 
calculation. Every point in z-direction is plotted and every other point in x- and j-direction is plotted. The 
simulation domain was bounded within a 42km x 42km x 18km domain. Mesh resolution was set to Ax = Ay = 
A z = 80/w. Time = Ylmin into the simulation. 
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Figure 49. The TASS computed Radar Reflectivity Factor field in the x-y plane at an altitude of 11120 m is 
shown in the left panel and the a An field on the same plane is shown in the right panel. Every other point in x- 
and j-direction is plotted. A Boeing 737-300 aircraft with a weight of \HHklbs is assumed in the a An 
calculation. Time = 37 min into the simulation. 
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VI. Summary 

Simulations of five severe convective induced turbulence events were performed using TASS. The events 
ranged from a large mesoscale convective system to isolated severe thunderstorms. The initial environments for 
TASS simulations were obtained from WRF predictions. Convection in TASS was initialized by an artificial 
impulse. A comparison of the model results with NEXRAD and other surface observations showed that it was able 
to simulate the evolution of storms in detail and with a good degree of accuracy. In most cases the model computed 
storm motion, storm orientation, peak echo tops, size of the radar echo, and the duration of the event agreed well 
with the observations but the peak radar reflectivity was generally over estimated. Some of the discrepancy in radar 
reflectivity is expected, since the model does not take into account the radar beam size and geometry in RRF 
calculation. 

Aircraft-centric turbulence characteristics were analyzed for each case. The turbulence analysis was consistent 
with previous studies (Proctor et al. 2002; Hamilton and Proctor 2002; Hamilton and Proctor 2003) and showed that 
large pockets of significant turbulence hazard can be found in regions of low radar reflectivity. Moderate and severe 
turbulence often was found in building cumulus turrets and overshooting tops, especially on the upwind side of the 
convective system. Data sets such as these can be very useful in radar simulation studies (e.g., Hamilton and Proctor 
2006) for testing and developing algorithms for hazard avoidance, as well as in the understanding of turbulence in 
relation to convective storms. 


Appendix A 

The Weather Surveillance Radar - 1988 Doppler (WSR-88D) is more commonly known as the NEXt Generation 
RADar or NEXRAD. The National Weather Service (NWS) operates the NEXRAD sites that cover the entire 
continental United States, Hawaii, Puerto Rico and Guam (http://www.roc.noaa.gov/WSR88D/). The NEXRAD 
data is archived by NCDC and is available for download from their website (http://www.ncdc.noaa.gov/nexradinv/). 
The NEXRAD base reflectivity data has a resolution of 1 km in range and 1° in azimuth. The data frequency in time 
is every 6min. The Level 2 (L2) data has 256 dBZ levels and the Level 3 (L3) data has 16 dBZ levels. The L2 base 
velocity is at a resolution of 250m x 1° compared to L3 base velocity which is at 1 km x 1°. The number of radar tilt 
angles is also larger in the L2 data (14 angles compared to 4 angles in L3). 

The NEXRAD Level 2 data is comprised of reflectivity, mean radial velocity and spectrum width, whereas the 
Level 3 data has 41 products which include base reflectivity, base velocity, base spectrum width, composite 
reflectivity at different levels, vertical integrated liquid, echo tops and wind profiles. The Level 3 products also 
include precipitation estimates such as one hour cumulative precipitation and total storm precipitation, etc. In 
addition, the Level 3 data contains derived information on storm structure and storm tracks, etc. The NEXRAD 
Level 2 and Level 3 products were used in the analysis and validation of TASS simulations presented in this paper. 
A brief description of the NEXRAD Level 3 products which were used in the analysis is given in Table A-l for 
reference. 


Table A-l: Description of Storm Characteristics Variables 


Variable 

Description 

Range 

Distance from radar in nautical miles. Maximum range = 124 nmi (nmi = 1.852km) 

Azimuth 

Azimuth angle from radar in degrees (0=north; 90=east; 180=south; 270=west) 

Echo Top 

Elevation of storm in kft (ft = 0.3048m) 

VIL 

Cell based vertically integrated liquid in kg/m 2 

Max Reflectivity 

Maximum reflectivity of the storm cell in dBZ from 0.5° elevation angle scan 

Storm Movement 

Storm movement in degrees and kts (kt = 0.5 14 m/s) 

One Hour Precipitation 

One hour precipitation in inches over the storm location (inch = 25.4mm) 
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